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Actions for dynamical fermion simulations: are we ready to go? 

K. Jansen'^* 

'^NIC/DESY Zcuthcn, Platanenallee 6, 15738 Zeuthen, Germany 

A critical review, playing devil's advocate, on present dynamical fermion simulations is given. 



1. Introduction 

No doubt, the lattice approach to quantum 
field theory has established itself as an integral 
part of high energy theoretical physics and has 
produced many results that are relevant to under- 
stand or even construct quantum field theories, 
as well as for phenomenology helping to interpret 
experimentally obtained data. Still, we are by 
far not at the end of the story. Serious dynam- 
ical fermion simulations with realistic number of 
quarks, i.e. two light and one heavy flavour, have 
just started and reveal already a number of unex- 
pected problems. It is the objective of this presen- 
tation to address these problems and discuss their 
relevance and possible consequences. Neverthe- 
less, we have made a lot of progress and we have 
understood much better to keep under control 
systematic errors such as discretization and finite 
size effects, the chiral limit and non-perturbative 
renormalization. However, still more patience has 
to be invested in our lattice computations to tell 
our colleagues from experiment trustworthy num- 
bers. New computer architectures are essential 
and will help tremendously to step forward. Of 
course, we are all hoping for a miraculous inven- 
tion of a new algorithm that would solve all our 
problems at once, although reality tells us that 
we will proceed in only small steps. 

2. Conceptual questions 

Certainly, for a non-lattice person (and even 
for a lattice practitioner) the plethora of actions 
that are used in todays lattice simulations is, to 
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say it moderately, confusing. The initial Wil- 
son fermion action T got many friends such as 
Symanzik improved actions 2,3^4 employing var- 
ious gauge actions such as the Wilson p], var- 
ious forms of Symanzik improved gauge actions 
|S] and renormalization group inspired actions 
A recent development concerning Wilson 
fermions is the formulation of twisted mass QCD 
j7j. Of course, there is the approach of staggered 
fermions [5] which comes again in improved forms 
"5* and with different gauge actions. Newer devel- 
opments are chiral invariant lattice fermions such 
as domain wall ^U], overlap JI] and perfect ac- 
tions 12 . If exact chiral symmetry is not insisted 
on, but only chirally improved fermions are con- 
sidered, designer actions such as FLIC ^SIj Hy- 
percube !14' or truncated fixed point fermions jl2j 
enter the stage. 

For me, the criterion for choosing an action 
should be guided by the demand that 

• the results obtained with these actions agree 
among themselves in the continuum limit 
and thus provide valuable cross checks 

• and if using a particular action will not 
waste computer time. 

There are a number of advantages and disad- 
vantages of the above mentioned actions, con- 
cerning field theoretical conceptual aspects, lat- 
tice artefacts, chiral and renormalization proper- 
ties and simulation cost that will be addressed 
below. Let me start by discussing locality prop- 
erties of the actions used. 

2.1. Rigorous actions 

The list of the above actions can be ordered 
by the level of rigor they obey in constructing a 
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lattice quantum field theory with the emphasis 
on their locahty properties. The initial, standard 
Wilson actions (including the twisted mass for- 
mulation without clover term) have been shown 
to be reflection positive jTJI, they have a posi- 
tive definite transfer matrix |16| and thus obey 
the reconstruction theorem (see first two refs. in 
|15| ) guaranteeing thus the reconstruction of the 
Minkowski Green's functions. Naive staggered 
fermions are on the same level |17) with the pe- 
culiarity that the transfer matrix has to be con- 
structed over two lattice spacings. 

2.2. Non-rigorous but local actions 

The next class of actions are the ones that are 
ultra-local, i.e. they have an interaction range 
of only a few lattice points, but for which reflec- 
tion positivity cannot be proven. To this class 
belong the Symanzik improved and the designer 
actions mentioned above. However, there cannot 
be a serious doubt that these actions are theoret- 
ically sound as the improvement terms added are 
irrelevant and vanish in the continuum limit. We 
might speak of "physical positivity" . 

2.3. Exponentially localized actions 

The chirally invariant lattice actions even give 
up ultra-locality jl8ll9| and only show an expo- 
nential localization |2()j , as shown for domain wall 
PT] . overlap [201 and perfect fermions [22] ■ Again, 
we have no reason to doubt that these actions 
define a local quantum field theory and describe 
QCD in the continuum limit. 

2.4. Non-local actions (?) 

Every time a lattice Dirac operator is con- 
structed, where the action ranges over all points 
on the lattice, the question of locality has to be 
asked anew. An example for such an operator 
are improved staggered fermions for exact, say, 
two flavour simulations ST. It is not clear, how 
such a staggered fermion operator can be defined 
and also in this contribution no solution will be 
provided. The suggestion that is used in exact dy- 
namical simulations, is to take the square root of 
the staggered operator [2SI , implemented by some 
polynomial approximation, which clearly defines 
an action that spreads over all lattice points. It 
should hence be tested, whether this operator is 



at least exponentially localized (see my plenary 
talk M for such a test on a Wilson operator) . 

There is a clear danger that through the square 
root operation a non- locality is introduced^. 

For a fixed set of parameters, i.e. fixed pion 
mass and fixed coupling, it should be expected 
that the localization of the operator is exponen- 
tial, with couplings decaying like e"''/'"'"'"'' with 
^'locai the localization range. In the continuum 
limit, staying at a fixed value of tq • m^, in prin- 
ciple two outcomes are possible. In scenario I, 
Hocai • "tItt = constant for a — > . In this case 
we would obtain a continuum field theory with 
a localization range of the order of the Compton 
wavelength of the pion. This would be clearly un- 
acceptable. In order to save such a theory at least 
to some extent, it would be necessary that the lo- 
calization range is much smaller than the Comp- 
ton wavelength of the heaviest physical particle 
of the target continuum field theory. 

In scenario II, riocai • — > for a . In 
this case we would get a local continuum quantum 
field theory. One may suggest a "wishlist" , i.e. to 

perform a check of the localization prop- 
erties of the square root staggered fermion 
operator. 

Such a work is in progress jZS.j . We assume first, 
by lack of an alternative, that the two fiavour 
staggered operator Mstaggercd is defined by tak- 
ing the square root of the four fiavor staggered 
operator. Second, we think of an exact simula- 
tion which will employ this square root opera- 
tor. In particular, following the arguments for 
deriving a strict locality bound of ref. 20 , we ex- 
pect the localization range to satisfy the bound 
^'locai • JTiquark = coustant. This bchaviour follows 
directly from considering admissible gauge field 
configurations, the knowledge of the maximum 
eigenvalue of the staggered operator and the fact 
that the smallest eigenvalue is determined by the 
square of the bare quark mass. It remains to be 



There might be additional shortcomings of such a defi- 
nition of a two flavour theory such as the violation of the 
optical theorem. However, it is not clear (at least not 
to me), how such a violation can be detected in practical 
simulations. 
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seen, whether this theoretical bound is satisfied in 
practical simulations or whether it grossly over- 
estimates the localization range as it was found 
in the overlap case |2U| . 

2.5. Perturbation theory 

For most of the actions just discussed, despite 
their sometimes complicated nature, perturba- 
tion theory has been worked out. There are fa- 
mous seminal papers for staggered fermions, see 
the first ref. in ^7] and as well as for Wil- 
son fermions '2T . These works build an impor- 
tant cornerstone to understand many aspects of 
lattice field theory such as renormalization. For 
more information on perturbation theory I refer 
to the recent review and the talk by H. Trot- 
tier at this conference |29| . 

The only aspect of perturbation theory I would 
like to mention here is the Reisz power counting 
theorem ,30 . This theorem guarantees that lat- 
tice integrals that appear in lattice perturbation 
theory exist in the continuum limit. The theorem 
is satisfied for Wilson fermions. However, stag- 
gered fermions do actually not obey one of the 
conditions of the theorem (see This does 

certainly not imply that there is a problem for 
staggered fermions in perturbation theory. How- 
ever, for a consistent and well founded perturba- 
tive discussion of staggered fermion it would be 
very helpful to 



construct a "Reisz theorem" for staggered 
fermions. 



Such a theorem for staggered fermions might 
also be very useful to analyze and underpin im- 
proved staggered actions, possibly leading to a 
disentanglement of the tastes in perturbation the- 
ory. 



The problem with our dynamical fermion sim- 
ulations is the "cost wall" . We know to some ex- 
tent the scaling behaviour of the algorithms and 
one formula for the cost C in the case of Wilson 
fermions according to ref. ,37, is 



C = A- 



(1) 



Here the exponents are given by z^r = 6, zl = 5 
and Za = 7; the value of the constant A can be 
found in |37| . Similar formulae are found by other 
groups Although this cost formula might 

still be plagued by a rather large uncertainty, it 
provides a good measure for the costs of the al- 
gorithms employed. 




Figure 1. A comparison of the cost of dynami- 
cal fermion simulations with (inexact) staggered 
(dotted line) and improved Wilson fermions (full 
line). The dashed line is the cost for Wilson 
fermions if the algorithms would perform a factor 
of four better than found in "37] . The left plot is 
for a value of the lattice spacing a = 0.09fm and 
the right plot for a — 0.045fm. 



3. Algorithmic questions 

Algorithms are the backbone of simulations in 
lattice field theory and testing and improving the 
algorithms cannot be appreciated enough. A re- 
view of existing algorithms for dynamical fermion 
simulations can be found in |32| . Promising new 
ideas can be found in ESI and see also 



In fig. n] we show a comparison of simula- 
tions with staggered (inexact, i.e. using the R- 
algorithm ^H]) and Wilson fermions. Several 
things are remarkable in this graph. The first 
is that for both kind of fermions at some value of 
the quark mass the costs explode. This is the fa- 
mous cost wall where the simulations become es- 
sentially un-feasible. The second point is that for 
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staggered fermions this point is reached at much 
lower values of the quark mass than for Wilson 
fermions. For the left plot we have used mea- 
sured cost data from improved staggered fermion 
simulations 39 and the cost formula eq. Q for 
Wilson fermions. A third observation is that the 
functional dependence of the cost on the quark 
mass seems to be the same for both actions, at 
least, given the accuracy of the data used here. 
The right plot shows again a cost comparison at 
a lower value of the lattice spacing, exhibiting a 
frightening need of teraflop years. The apparent 
advantage of the simulation cost using improved 
staggered fermions makes it even more important 
and urgent to test the theoretical basis of this 
approach as discussed in section 2. 

One shortcoming of the staggered fermion sim- 
ulations so far is that only a non-exact algorithm 
is used 10 . In principle we know how to ex- 
ecute an exact odd flavour algorithm A 
comparison of the in-exact algorithm and the ex- 
act one is shown in fig. El A striking effect of 
this graph is that the extrapolation to zero step 
size as needed in the in-exact simulation is non- 
monotonic. Hence such an extrapolation has to 
be done with very great care using many simu- 
lation points. Since exact odd flavour algorithms 
are on the market, it should be urged that they 
are indeed used in the staggered fermion simula- 
tions. 

As a conclusion, it would be very useful to 
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use and test exact algorithms 
and to perform a fair cost comparison of 
exact algorithms including the approach to 
the continuum limit. 



The actions discussed in section 2 are some- 
times very complicated, including often fatten- 
ing of links. A good news is that methods have 
been developed that allow for simulations of such 
actions either by reject /accept steps or by 
computing the force explicitly 02| using special 
projection methods. 

As a last topic in this algorithm section I would 
like to discuss a problem of principle. Random 
matrix theory (RMT) 03] predicts the eigenvalue 
distribution for the lowest non-vanishing eigen- 
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0.00 0.02 0.04 0.06 0.08 0.10 
dt 

Figure 2. A comparison of the plaquette from an 
exact odd flavour simulation (triangles) and an 
in-exact one (circles, square) against the step-size 
dt, from the first reference of |41j . 



value A in topological charge sector zero to be 
P(A) = |e-3^' with z = AST/ and S the infinite 
volume scalar condensate This prediction is 
confirmed in practical simulations j45| . 

Such a distribution predicts the occurrence of 
very small eigenvalues to be frequent, and hence 
close to the critical point huge fluctuations have 
to be expected, rendering the simulations very ex- 
pensive. Of course, a non- vanishing quark mass 
will regulate the problem, but it can do so only 
to a limited amount when the physical point at 
realistic values of the quark mass is reached. The 
conclusion is therefore that we have to expect 
that dynamical fermion simulations are becoming 
even more problematic than the cost formula of 
eq. would imply when approaching small quark 
masses. As a result, the need to make contact 
with chiral perturbation theory is most urgent. 

4. Scaling in the quenched approximation 

Before turning to any results of unquenched 
simulations, let me discuss scaling in the quenched 
case first. Surprising results were presented by 
Aoki 46 at Lattice2000. He used different ac- 
tions: staggered fermions with plaquette action, 
Wilson and improved Wilson fermions with pla- 
quette and renormalization group inspired gauge 
actions. His results for mjq/mv for the different 
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actions as a function of mpg/my (i.e., the Edin- 
burgh plot) did not agree in the continuum limit 
even for for heavy quark masses. The staggered 
results disagreed with the Wilson results, which 
agreed among themselves. Surprisingly, there was 
not much response to this striking and worrisome 
observation afterwards. 

The question arises, whether something could 
be wrong with the continuum limit of improved 
staggered fermions. In Aoki's article, a partic- 
ular form of the continuum extrapolation was 
taken, linear in the lattice spacing for Wilson type 
fermions and an dependence for the staggered 
results. In fig. Ell show an alternative fitting pro- 
cedure, allowing for an additional a^-dependence 
for the Wilson plaquette data'^. 
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Figure 3. Continuum extrapolations of unim- 
proved Wilson and unimproved staggered 
fermions, employing a pure dependence in the 
unimproved staggered and an a -I- dependence 
in the unimproved Wilson case. A fixed value 
of mps/my = 0.7 was chosen. Also included 
are data for improved staggered, hypercube and 
FLIC fermions. The data were compiled together 
with J. Zanotti. 



Although this way of fitting the data might not 
be the final truth, it provides a possibility with 

leave out the improved Wilson data since they arc not 
very consistent among themselves and it would certainly 
be important to generate a new set of data. 



the conclusion that there is nothing wrong with 
the continuum limit of staggered fermions. The 
continuum value for mpf/my is also consistent 
with results for various other improved actions as 
included in the plot. It would be very important, 
to supplement more and better data for such a 
plot and to 



perform a precise scaling analysis for var- 
ious fermion actions in the quenched ap- 
proximation. 



Such a test would be very re-assuring that all 
the calculations that are done with the many dif- 
ferent actions lead to consistent values in the con- 
tinuum limit and can hence lead to a good control 
over the systematic errors. For attempts for such 
a precise scaling tests I refer to the figure (done 
together with J. Zanotti) in my plenary talk [23 
and the talks by A. and P. Hasenfratz at the con- 
ference EaiH]. 

5. Problems with dynamical fermion sim- 
ulations 

I now turn to somewhat surprising problems 
that are already encountered in practical simu- 
lations with dynamical fermions performed as of 
today. In the light of these problems, it may not 
be wise to start unthinkingly large scale and very 
expensive simulations with dynamical fermions 
but rather to concentrate on these problems for a 
while and try to understand and solve them. 

5.1. Gauge actions 

As mentioned before, improved gauge actions 
are not reflection positive. This leads to complex 
eigenvalues in the transfer matrix 5 . In fig.^the 
static potential V{ro) is shown at a distance of the 
hadronic scale tq ^HI as function of the euclidian 
time. The values of V(ro) for the improved gauge 
actions approach their (in time) asymptotic value 
from below, indicating the effects of the complex 
eigenvalues [SO]. Although for asymptotic times 
the value of the potential can be determined, this 
suggests that for these actions larger time extents 
would have to be used. This can have conse- 
quences in particular for determinations of glue 
ball masses from such actions. 
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mal cycles^, clear indications of a first order phase 
transitions were detected as can be seen in fig. 
where a pronounced hysteresis effect shows up. 
A similar effect appears also for Nf = 2 flavours 
|54| if again the plaquette gauge action and non- 
perturbatively improved Wilson fermions 0] are 
used. It might be speculated that these phase 
transitions are related to the phase structure 
of the fundamental-adjoint pure gauge action. 
Working in the vicinity of such a phase transi- 
tion can lead, e.g., to unexpectedly large lattice 
artefacts that may render the continuum extrap- 
olation difficult. 



Figure 4. Static potential for various gauge ac- 
tions. Note that, as the effect of complex eigen- 
values, in the case of improved gauge actions [H] 
the potential starts from below. 



From a free field analysis (50) it also appears 
that different gauge actions might lead to very 
different lattice artefacts. In addition, the lattice 
spacing dependence need not to be very smooth 
approaching the continuum limit thus driving 
a reliable continuum extrapolation difficult. It 
might be useful therefore to perform simulations 
with different actions to check the continuum 
limit. An example, from my own experience, for 
such a procedure can be found in ref. [IJ where 
simulations using different fermion actions were 
performed to extract moments of parton distribu- 
tion functions. We will come back to this point 
towards the end of this proceedings. 

Another, last difficulty of improved gauge ac- 
tions is that they are very inefficient to sample 
topological charge sectors Although this 

is not too problematic for quenched simulations 
since there many configurations can be generated 
fast, it may have serious consequences for dynam- 
ical fermions when such actions are employed. 

5.2. Wilson fermions 

The CP-PACS collaboration started Nf = 3 
simulations with improved Wilson fermions and 
the plaquette gauge action. By performing ther- 
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Figure 5. Thermal cycles in k for the plaquette 
for Nf = 3 non-perturbatively improved Wilson 
fermions at various values of f3 ranging from P = 
4.6 to P — 6 from bottom to top (from [5^1 )■ For 
f3 5 hysteresis effects are observed, indicating 
the existence of a first order phase transition. 



Using the Iwasaki action [H' or the Symanzik 
gauge action the signs of the first order phase 
transitions seem to disappear completely |53| . 
However, it is not clear, whether such effects even- 
tually will strike back. Therefore, it would be 



*In this standard procedure to detect phase transitions, 
one starts at some value of Kstart, increases the value of 
K by some Aft until a value of ftgnai is reached and then 
runs back to Kgtart- 
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important to 



study and understand the T — phase 
diagram of QCD, to investigate different 
gauge actions and to understand the na- 
ture of the phase transition. 




Figure 6. The quark mass difference of m(oo) — 
m(L) in quenched and in full QCD using non- 
perturbatively improved Wilson fermions [L/a = 
oo = L / a = 16). Note that for the unquenched 
case the lattice artefacts are substantially larger. 
As an aside, we also show results from standard 
Wilson fermions (W-W) where the lattice arte- 
facts are even bigger (from |55]). 



In fig. ini we show a comparison of the lat- 
tice artefacts in the quark mass for quenched 
and unquenched simulations j^S] using non- 
perturbatively improved Wilson fermions. The 
figure indicates that for dynamical fermions the 
artefacts are substantially larger, although they 
are still smaller than using naive Wilson fermions 
|56j . Whether this increase of lattice artefacts 
is due to the possible phase transition discussed 
above or originates from a completely different 
source, is an open question and gives additional 
motivation to understand the phase diagram bet- 
ter. 

During the conference it has been pointed out 
|57| that in dynamical simulations using Wil- 



son or staggered fermions the hadronic scale rg 
depends rather strongly on the quark mass at 
fixed values of /3. Clearly, this effect has to be 
taken into account when quantities are analyzed 
as functions of, say, rp • mps- This becomes es- 
pecially important when comparisons with chiral 
perturbation are performed. In addition, it was 
shown that in the quantity ro(m)/ro(mrcf) with 
TTircf determined by mps/my — 0.45, staggered 
fermions exhibit surprisingly large lattice arte- 
facts. Clearly, it would be important to check 
and clarify this finding further with results from 
improved staggered fermion simulations. 

5.3. Staggered fermions 

For simulations with staggered fermions, naive 
or improved, there is a problem with the topo- 
logical charge identification when coarse lattices 
at large values of the coupling are used. This 
can be seen quite evidently in studies that com- 
pare the predictions of the eigenvalue distribu- 
tions predicted by RMT with the results of nu- 
merical staggered fermion simulations. Data ob- 
tained in topological charge sector 1 are fitted 
with the topology prediction of RMT, while the 
curves for topology 1 fail completely These 
studies were mainly done with unimproved stag- 
gered fermions. For newer results with improved 
staggered fermions, see |47j . 

In 1^ it was argued that this failure of stag- 
gered fermions with respect to questions related 
to topology has to occur when the gauge fields 
are not smooth enough and one is not working at 
small enough values of the coupling. The reason is 
that with such rough gauge fields, the taste break- 
ing effects are strong and the symmetry breaking 
pattern is not the one of the continuum. 

Indeed, in the first reference of '58j where sim- 
ulations of the Schwinger model have been per- 
formed, there are indications that for larger val- 
ues of f} the predictions of RMT and the numeri- 
cal data get reconciled. Nevertheless, it would be 
very important 



to check explicitly that the predictions of 
RMT and other properties of topology are 
reproduced by, preferably, improved stag- 
gered fermions. 
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6. Results for dynamical Wilson fermions 

This section is devoted to some results that 
have been obtained recently with dynamical 
Wilson fermions. Results concerning staggered 
fermions can be found in |60 | and I will only 
shortly discuss staggered fermions in the context 
of chiral perturbation theory here. 

6.1. Chiral perturbation theory 

The algorithms we are using at the moment will 
sooner or later hit a cost wall where the exponen- 
tial growth of the cost will set in with full power 
when the quark mass is lowered below a threshold 
value. In order to obtain results for physical val- 
ues of the pion and p masses it becomes therefore 
mandatory to make contact to chiral perturba- 
tion theory (xPT) and then use xPT to perform 
the extrapolations to the "physical point" . Find- 
ing out the window where on the one hand lattice 
simulations can be performed and, on the other 
hand, xPT can be controlled and higher order 
effects can be safely neglected, is of central im- 
portance in lattice QCD and much effort will be 
invested to resolve this question. Fig. 1 suggests 
that improved staggered fermions may reach this 
overlap region sooner than Wilson fermions. 

There are two main strategies to confront lat- 
tice results with predictions of xPT: 

• the first strategy is to extrapolate lattice 
data to the continuum limit first and then 
compare to xPT; 

• the second strategy is to take the lattice 
artefacts into account in xPT itself and 
compare then to lattice data. 

Let me discuss the first strategy first. The big 
advantage of this procedure is that in the con- 
tinuum chiral symmetry is restored and therefore 
a direct comparison to xPT becomes possible^. 
Clearly, the bad side of the coin is that perform- 
ing the continuum limit is demanding and needs 
large computer resources. 

This computational demand is at present 
"tricked out" by taking lattice data at only one 

^Of course, using chiral invariant lattice fermions would 
also fulfill this requirement, but, they are, unfortunately, 
very expensive to simulate. 



value of the lattice spacing. Of course, the im- 
proved lattice actions in both the Wilson and the 
staggered cases substantially reduce the lattice 
artefacts. However, only future simulations can 
tell, whether the lattice data were indeed close 
enough to their continuum values. 

In ref. (BTj a recent analysis of chiral pertur- 
bation theory applied to lattice data are given. 
Older discussions are numerous and I have to re- 
fer to "W?, the plenary discussion at Lattice 2002 
and 63 for other examples (taken from the anal- 
ysis of moments of parton distribution functions, 
my own hobby-horse). 
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Figure 7. Comparison of CP-PACS data with dif- 
ferent kinds of performing xPT (from The 
dashed dotted curve is a parameter-free third or- 
der result. The dotted line is a fourth order rep- 
resentation using a best fit for the additional low 
energy parameters appearing in the fourth order 
of xPT. 



In the formulae used to perform the analyzes 
in [ST], the parameters of the chiral Lagrangian 
were fixed by taking experimental input values 
obtained at the physical point. In the calcula- 
tion chiral symmetry was consistently kept and 
an "improvement term" to cancel cut-off effects 
was added. The numerical data are taken from 
improved Wilson CP-PACS simulations and are 
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only at one value of the lattice spacing. Fig. [7| 
shows a comparison of the predictions from xPT 
in this setup and the data, which seems to be not 
so bad, even for quite large pion masses. However, 
in ref. ^61 itself it is claimed: We stress again that 
applying the expressions to pion masses above 600 
MeV is only done for illustrative purposes, for a 
realistic chiral extrapolation smaller pion masses 
are mandatory. 

The reason for this becomes clear by looking 
at the size of corrections that can appear in a 4- 
loop calculation. These corrections can be very 
large, depending on the precise form of perform- 
ing xPT. I do not have the space (and it is not the 
goal of this article) to discuss in detail which kind 
of xPT is correct, what flaws there might be and 
where conceptual difficulties appear. However, I 
believe, that finally it has to be pure xPT that 
has to be used, without any additional modeling 
or assumptions. Only then can we safely employ 
xPT to do for us the job of extrapolating to the 
physical point. 

The second strategy is to take possible lattice 
artefacts into account in xPT itself. This idea 
was put to the fore by S. Sharpe 1641 and was 
then further developed in (>•") for Wilson and in 
[HH| for staggered fermions. The advantage of this 
approach clearly is that chiral symmetry violat- 
ing effects, taste violating effects in the case of 
staggered fermions and even different actions can 
be taken into account from the beginning. The 
disadvantage is that the corresponding chiral La- 
grangian has twice as many low energy constants, 
and that the additional parameters depend on the 
bare coupling go leading to a plethora of fit pa- 
rameters. 

This approach has been used already in practi- 
cal simulations. For results with improved stag- 
gered fermions I refer to the contribution of S. 
Gottlieb at this conference (^01. The only re- 
mark, I want to add here is that xPT provides 
a justification of taking the square root of the 
fermion determinant. This situation corresponds 
to a partially quenched setup in which two quarks 
are taken to be quenched j^. The idea can be 
generalized using the replica trick: take arbitrary 
numbers Nu, Nd, Ns for the u, d and s quarks. 
Perform the computation in this general form and 



set in the end ^ Nj, ^ Ns = 1/4 [HHI- In 
although in a different context, it has been 
demonstrated that such a procedure is, at least to 
1-loop order, completely equivalent to the super 
symmetric method and, although a proof to all 
orders is still missing, there is a good chance that 
the argument might work. 

After this excursion to staggered fermions, let 
me now turn back to Wilson fermions. There are 
two works at the moment where simulation data 
are confronted to the predictions of xPT when 
lattice artefacts are included in xPT. The first 
work is ref. [201 and the second is ref. ^V^. We 
show in fig. ISI a plot of mpg/mq as a function 
of the quark mass at various values of (3. For 
each value of /3 two curves are shown. The curves 
showing a very strong curvature represent con- 
tinuum xPT while the lines that fit the data rep- 
resent xPT with lattice artefacts taken into ac- 
count. Clearly, continuum xPT is not able to de- 
scribe the simulation data while lattice xPT fits 
the data well. Another observation is, however, 
that very large lattice artefacts in rripg/mq seem 
to exist. Data at various values of f3 come out to 
be quite different. 

Another comparison is shown in fig. |3 (from 
\7l\). Here a ratio of decay constants of the type, 
fvs/ifwfss) is plotted with V and S denot- 
ing valence and sea quantities. In this double 
ratio, which is only one example of the many fits 
that can be found in ref. [7T], a good agreement 
with the predictions of xPT is seen. However, 
there seem to be tremendous cancellations of cut- 
off effects, keeping in mind that the simulations 
have been done at a value of the lattice spacing 
of a = 0.28fm. 

Another, somewhat worrisome observation is 
that the universal A-parameters come out to be 
very different, depending on the value of /3 that 
is used (in case of jZOI) or how the lattice spacing 
is defined (in case of |7J). Clearly, we are just in 
the beginning of using the new approach of xPT 
where lattice spacing effects are taken into ac- 
count and the first results, using this method, are 
somewhat puzzling. In particular, I would like to 
advocate to de-double the double ratios in order 
to be able to disentangle the cut-off effects and 
show their size, furthering thus our understand- 



10 




0.05 0.1 0.15 

^Awi (GeV) 



0.2 0.25 




0.4 o.e o.f 



Figure 8. rripg/mq as a function of the quark 
mass for various values of /3. For each value of /3 
two curves are shown. The curves showing a very 
strong curvature represent continuum xPT while 
the lines that fit the data represent xPT with 
lattice artefacts taken into account, from [70]. 



ing of xPT and the connection to lattice QCD. 
7. Numerical results 

In this section, I would like to discuss a few 
selected results from recent dynamical Wilson 
fermion simulations. 

7.1. Finite size effects 

It seems that we finally see the exponential 
finite size effects of hadron masses. The gen- 
eral form of these finite size effects are M(L) — 

M = -Ta;^r^^^(^^/)eV^^ where F is 
the TT — TT forward scattering amplitude in in- 
finite volume |75j- A calculation of the inte- 
gral leads to the exponential finite size correc- 
tion M{L) - M (X L-3/2g-m„i^ Corrections to 
this formula have been computed within xPT 
[73] ■ An alternative calculation was performed 
and also tested in ref. [73]. At this conference 
there were two contributions |74I75| that showed 
that indeed the finite size effects are exponential 
and not power like of the form M = rrioo + c/L^ 
as anticipated in ref. [TSJ. In fig. ^] a compar- 
ison of numerical simulation data to various fits 



Figure 9. The double ratio mentioned in the text 
compared to xPT. ^ denotes the ratio of valence 
to sea quark mass. 



of the finite size effects are shown. It seems that 
indeed the exponential form of finite size correc- 
tions are preferred by the data confirming thus 
the exponential finite size effects in lattice simu- 
lations. A similar picture and conclusion can be 
found in "TC. Clearly, these first results for the 
finite size effects need to be confirmed in future 
simulations. 

7.2. Fundamental parameters of QCD 

The running coupling and running quark 
masses are basic and important parameters of 
QCD that can be extracted from lattice simula- 
tions. The status of these calculations were re- 
ported at this conference [ZJ]. Although the run- 
ning itself has been determined already to a quite 
high precision, the determination of the physical 
scale is still missing. In addition, it would be 
necessary to complement the present simulations 
with even larger lattices to keep the continuum 
limit under better control. Both of these goals 
will certainly be done in the near future. In the 
past, the Schrodinger functional, which was used 
for the scale dependent renormalization [JH], was 
the domain of the ALPHA collaboration. Nowa- 
days, there are new groups who use this renormal- 
ization scheme. In particular, a serious attempt 
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Figure 10. Finite size effects of the pion, p and 
Nucleon masses compared to various fits corre- 
sponding to different analytical predictions (from 
|75|). The exponential form of finite size effects 
seem to be clearly preferred. 



for Nf — 3 flavours of dynamical fermions has 
been started by the CP-PACS and JLQCD col- 
laborations 179'. 

The CP-PACS and JLQCD collaborations 
seem not only to join their forces to compute 
the coupling and quark masses but also for other 
quantities. This is a convincing indication that 
dynamical fermions are a real challenge and a 
very difficult problem. Present simulations em- 
ploy the Iwasaki gauge action of ref. W combined 
with the Symanzik improved Wilson fermion ac- 
tion. The reason for this choice is the occurrence 
of a phase transition when the standard gauge ac- 
tion is employed as mentioned above. The value 
of csw is the non-perturbatively improved one 

In the Nf = 3 simulations, the physical extent 
of the lattice is about L « 1.6 fm. With such a 
lattice, it is not possible to perform simulations 
for baryonic quantities reliably since the finite size 
effects will be very large Hence, at the mo- 
ment one concentrates on mesons. In fig. 1111 the 
relative difference of the K and $ meson masses as 
compared to the experimental data is shown. The 
open triangles represent the quenched approxi- 
mation where for different physical input value a 



large variation is found. It is very reassuring that 
for dynamical fermions this intrinsic uncertainty 
of the quenched approximation seems to be elim- 
inated completely. 



m^. (K-input) (K-input) ((|)-input) m^. {(^-input) 



_ «N,.3 




i * 




■ i 

: 5 


— ; 

- 5 - 







Figure 11. Relative difference of the K and $ 
meson masses as compared to the experimental 
data. Open triangles represent the quenched ap- 
proximation. 



Similar to what has been said above for the 
quenched approximation, it would be very inter- 
esting to 

combine Wilson, staggered and other dy- 
namical fermion action results to check the 
different approaches in the chiral, contin- 
uum and infinite volume limits. 



7.3. Other topics 

An article like the present one has to be selec- 
tive and a number of topics had to be left out, un- 
fortunately. In the first place, I could not discuss 
progress with dynamical domain wall jg2'] or even 
overlap |83j fermions. For domain wall fermions it 
seems that an appropriate gauge action has been 
determined where an additional adjoint piece is 
added which seems to solve some problems with 
unexpected large values of the residual mass. 
The simulations for producing physics have been 
started and I am sure that at the next confer- 
ence we will see a lot of interesting results. I 
would also like to mention some works that dis- 
cuss the study of violations of CP |H1] and even 
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CPT |85l86j without being able to go into any 
detail. 

I could also not go into results for moments 
of parton distribution functions from dynamical 
simulations 87. , nor could I discuss questions 
concerning topology in depth jHS]. Attempts to 
tackle the notoriously difficult problem of the 77' 
meson by making use of the eigenvalue spectrum 
can be found in |89j . My personal opinion here 
is that this is an interesting approach, but it still 
has a number of open questions that should be 
addressed in the future. 

8. Conclusion 

Obviously, there are dangerous animals on the 
lattice 90 . Our present dynamical fermion sim- 
ulations are still in a rather early stage but 
they meet already such dangerous animals: un- 
expected phase transitions, questions of field the- 
oretically non-local actions, actions with unitarity 
violations, occurrence of large lattice artefacts, al- 
gorithms that become extremely costly when ap- 
proaching the physical point, uncertainties in the 
continuum limit of different lattice actions, diffi- 
culties with the chiral extrapolations. 

In the real world, we have learned to tame lions 
and tigers and in lattice field theory this will be 
done for sure, too. However, we have to spend 
work for achieving this: we have to explore and 
understand the phase diagram of (zero tempera- 
ture) QCD, we have to study the locality proper- 
ties of taking the square root, we have to explore 
better the phase diagrams of extended gauge ac- 
tions and understand the properties of such ac- 
tions. We should perform scaling tests first in 
the quenched approximation to demonstrate to 
ourselves and to the rest of the world that uni- 
versality holds (with all the caveats of the notion 
of universality in the quenched approximation). 
Participate in algorithm development. It is here, 
where a breakthrough can happen which would 
solve the cost- wall problem we have with existing 
algorithms; it would also be very helpful to join 
the ILDG [HI] initiative to combine our efforts 
and exchange configurations. 

One recommendation could be to not put all 
cards on one particular action. In fact, it would 



be wise to first select carefully actions that do 
not show the flaws and problems mentioned in 
this contribution. Even after such identifications 
have been made, it may be worthwhile to perform 
simulations with two (good) actions to check for 
systematic errors. 

Another recommendation could be to use ac- 
tions tailored for the physics problem one is in- 
terested in. If problems with rather heavy quark 
masses are considered, (theoretically sound) stag- 
gered fermions or improved Wilson fermions are a 
good choice. Going lighter, designer actions such 
as FLIC, hypercube, truncated perfect actions or 
domain wall fermions with rather small extra di- 
mension might enter the game. Finally, for really 
light quarks, chiral invariant (up to machine pre- 
cision) domain wall and overlap fermions should 
be used. 
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